Geometry of the quasi-hyperbolic Szekeres models 



Andrzej Krasinski 
N. Copernicus Astronomical Centre, 
Polish Academy of Sciences, 
Bartycka 18, 00 716 Warszawa, Poland 



(N 

o 

< 
ov 



Krzysztof Bolejko 
Sydney Institute for Astronomy, School of Physics A28, 
The University of Sydney, NSW 2006, AustraliJ^ 
(Dated:) 

Geometric properties of the quasi-hyperbolic Szekeres models are discussed and related to the 
quasi-spherical Szekeres models. Typical examples of shapes of various classes of 2-dimensional 
coordinate surfaces are shown in graphs; for the hyperbolically symmetric subcase and for the 
general quasi-hyperbolic case. An analysis of the mass function M(z) is carried out in parallel 
to an analogous analysis for the quasi-spherical models. This leads to the conclusion that M(z) 
determines the density of rest mass averaged over the whole space of constant time. 

PACS numbers: 
Keywords: 



> 

O 

(N 
00 

O 



X 



I. MOTIVATION 

Continuing the research started in Refs. [l[ and 
the geometry of the quasi-hyperbolic Szekeres models is 
investigated. Unlike the quasi-spherical Szekeres mod- 
els that have been extensively investigated Q - [2(| and 
are rather well-understood by now, the quasi-plane and 
quasi-hyperbolic models are still poorly explored. This 
situation has somewhat improved recently: in Ref. [l[ a 
preliminary investigation of the geometry of both these 
classes was carried out, and in Ref. 0] it was shown that 
the physical interpretation of the plane symmetric mod- 
els becomes clearer when a torus topology is assumed for 
the orbits of their symmetry. 

The present paper is an attempt to understand the ge- 
ometry of the quasi-hyperbolic model. In Sec. |TT] the 
full set of the j3, z ^ Szekeres solutions is presented. 
In Sec. IIIII limitations for the arbitrary functions in the 
quasi-hyperbolic models are discussed that result from 
the spacetime signature and from the evolution equa- 
tion. It is also shown that a set where the mass func- 
tion is zero is allowed to exist. In Sec. IIVI it is recalled 
after Ref. that the quasi-hyperbolic Szekeres mani- 
fold is all contained within an apparent horizon, i.e. is 
globally trapped. In Sec. IVII the geometry of various 
2-dimensional surfaces in the hyperbolically symmetric 
subcase is investigated and illustrated with graphs. In 
Sec. IVIII it is shown what deformations to the surfaces 
of constant t and if ensue in the general quasi-hyperbolic 
case. In Sees. IVIIII - IXII various properties of the mass 
function in the quasi-spherical models are discussed, in 
order to prepare the ground for an analogous discussion 



of the quasi-hyperbolic case. This latter task is carried 
out in Sees. IXIII and IXIIII The purpose of this was to 
identify the volume in a space of constant t, which could 
be related to the mass M{z). This goal was not achieved 
as intended, but it was shown that M(z) determines the 
density of rest mass averaged over the space of constant 
time. Section IXIVI is a summary of the results. 

The aim of this paper is to advance the insight into the 
geometry of this class of spacetimes. This is supposed to 
be the next step after the exploratory investigation done 
in Ref. [l|. 



II. INTRODUCING THE SZEKERES 
SOLUTIONS 

This section is mostly copied from Ref. |2j, mainly in 
order to define the notation. 

The metric of the Szekeres solutions is 



ds 2 = dt 2 



^ (dz 2 +dy 2 



(2.1) 



where a and (3 are functions of (t, x, y, z) to be deter- 
mined from the Einstein equations with a dust source. 
The coordinates of (|2.1[) are comoving, so the velocity 
field of the dust is = S^o, and = 0. 

There are in fact two families of Szekeres solutions, 
depending on whether (3, z = or f3, z ^ 0. The first family 
is a simultaneous gene ralisation of the Friedmann and 
Kantowski - Sachs [2l| models. Since so far it has found 
no useful application in astrophysical cosmology, we shall 
not discuss it here (see Ref. [17] ) , and we shall deal only 
with the second family. 

After the Einstein equations are solved, the metric 
functions in (12.11) become 
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c 73 = $(t, z)e v(z - x ' y) , 

e a = h(z)<S>(t,z)(3, z =h(z)(<S>, 



(2.2) 
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e~" = A(z)(x 2 +y 2 ) +2B 1 {z)x + 2B 2 {z)y + C(z), 
where $(t, z) is a solution of the equation 



-k(z) 



2M(z) 



-A$ 2 



(2.3) 



while h(z), k(z), M(z), A(z), B x (z), B 2 (z) and C(z) are 
arbitrary functions obeying 



g( Z ) = HAC- 



B X 2 ~B 2 2 ) = l/h 2 (z) + k(z). (2.4) 



The mass-density p is 



(2Me 3 ") , 2 



o2/3 f p /3 



K = 8ttG/c 4 



(2.5) 



This family of solutions has in general no symmetry, 
and acquires a 3-dimensional symmetry group with 2- 
dimensional orbits when A, B±, B 2 and C are constant 
(then v lZ — 0). The sign of g(z) determines the geometry 
of the surfaces of constant t and z and the symmetry of 
the v, z = subcase. The geometry is spherical, plane or 
hyperbolic when <7>0,<7 = 0org<0, respectively. 
With A, B\, B 2 and C being functions of z, the surfaces 
z = const within a single space t = const may have dif- 
ferent geometries, i.e. they can be spheres in one part 
of the space and surfaces of constant negative curvature 
elsewhere, the curvature being zero at the boundary - see 
a simple example of this situation in Rcf. [I]. 1 The sign 
of k(z) determines the type of evolution when A = 0; with 
k > the model expands away from an initial singularity 
and then recollapses to a final singularity, with k < the 
model is ever-expanding or ever-collapsing, depending on 
the initial conditions; k — is the intermediate case with 
expansion velocity tending to zero asymptotically. 

The Szekeres models are subdivided according to the 
sign of g(z) into quasi-spherical (with g > 0), quasi-plane 
(g = 0) and quasi- hyperbolic (g < 0). The geometry of 
the latter two classes has, until recently, not been in- 
vestigated and is not really understood; work on their 
interpretation was only begun by Hellaby and Krasihski 
[l|, and somewhat advanced for the quasi-plane models 
by the present author Q. The sign of g(z) imposes lim- 
itations on the sign of k(z). For the signature to be the 

physical (-1 ), the function h 2 must be non-negative 

(possibly zero at isolated points, but not in open subsets), 
which, via (|2.4p . means that g(z) ~ k(z) > everywhere. 
Thus, with g > all three possibilities for k are allowed; 
with .9 = only the two k < evolutions are admissible 
(k = only at isolated values of z), and with g < 0, only 
the k < evolution is allowed. 

The quasi-spherical models may be imagined as such 
generalisations of the Lemaiitre - Tolman (L-T) model 



1 In most of the literature, these models have been considered sep- 
arately, but this was only for purposes of systematic research. 



in which the spheres of constant mass are non-concentric. 
The functions A(z), B\(z) and B 2 (z) determine how the 
center of a sphere changes its position in a space t = 
const when the radius of the sphere is increased [l6| . 

Often, it is practical to reparametrise the arbitrary 
functions in the Szekeres metric as follows [22j ■ Even if 
A = initially, a transformation of the (x, y)-coordinates 
can restore A ^ 0, so we may assume with no loss 

of generality p7j |. Then let g ^ 0. Writing 



(A,B U B 2 ) 



25_ 

k = -\g\x2E, M 



l9l (l,-P,-Q), e d ^g/\g\, (2.6) 



\g\ 3/2 M, $ 

we can represent the metric ()2.1|) as 



\g\R, 



C V def „ def S 

~= ~ 2 



1.91 



x - P 
S 



ds 2 =dt 2 - {R > z - R£ ': /£) dz 2 ~ 

e + 2E(z) 



S 

E. 

E 2 



, (2-7) 



(da; 2 + dy 2 ) . 



(2.8) 



When g = 0. the transition from (|2~Tj) to (|2~7l> - (|2~8| is 

A= 1/(25), B l = -P/(2S), B 2 = -Q/(2S), k = -2E, 
M = M and $ = R. Then JH) - dH]) applies with 
e = 0, and the resulting model is quasi-plane. 
Equation (|2.3|) . in the variables of (|2.8p becomes 



R,t 2 = 2E(z) 



2M(z) 
R 



(2.9) 



From now on, we will use this representation. The for- 
mula for density in these variables is 



npc 



2 (M» -3M£, Z /£) 
R 2 (R, z -R£, z /£) ' 



(2.10) 



For p > 0, (M, z -3M£ :Z /£) and (R, z -R£, z /£) must 
have the same sign. Note that the sign of both these 
expressions may be flipped by the transformation z — > 
— z, so we may assume that 



R, z -R£ lz /£>0 



(2.11) 



at least somewhere. In this preliminary investigation we 
assume that we are in that part of the manifold, where 
(|2~TTD holds. 

In (|2.7p - (|2.8p the arbitrary functions are indepen- 
dent. 2 However, (|2.7p - (|2.8p creates the illusion that 
the values e = +1,0,-1 characterise the whole space- 
time, while in truth all three cases can occur in the same 
spacetime. 



Equation |g3) defines C = ylsl [(P 2 + Q 2 ) /S + eS] /2. 
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Within each single {t = const, z — const} surface, in 
the case e — +1, the (x,y) coordinates of (|2.ip can be 
transformed to the spherical ($, ip) coordinates by 

(x - P,y- Q)/S = cot(t?/2)(cos<^,sm<p). (2.12) 

This transformation is called a stereographic projection. 
For its geometric interpretation and for the correspond- 
ing formulae in the e < cases see Refs. Q and_[l7j|. 
The shear tensor for the Szekeres models is [l7( 



E = 



= -E diag (0,2,-1,-1) 

R,tz —R,t R,z I R 
R, z —R£,z /£ 



where 



(2.13) 



Since rotation and acceleration are zero, the limit o~ a p — > 
must be the Friedmann model [l?], HH . In this limit we 
have 



R(t,z)=r(z)S(t), 
and then (12.91) implies that 

E/r 2 ,M/r 3 and tsiz) are all constant. 



(2.14) 



(2.15) 



However, with P(z), Q(z) and S(z) still being arbitrary, 
the resulting coordinate representation of the Friedmann 
model is very untypical. The more usual coordinates re- 
sult when 



R',z Q?z ^s; 



o. 



(2.16) 



and r(z) is chosen as the z'-coordinate. (We stress that 
this is achieved simply by coordinate transformation, but 
writing it out explicitly is an impossible task.) However, 
(|2~T4]) and (|2~T^]) substituted in (|2.8p give the standard 
representation of the Friedmann model only when e — 
+1. With e = and e = — 1, further transformations are 
needed to obtain the familiar form [l|, [24|, HH] . 

The above is a minimal body of information about the 
Szekeres models needed to follow the remaining part of 
this paper. More extended presentations of physical and 
geometrica l prope rties of these models can be found in 
Refs. Q, II [Mil- 



III. SPECIFIC PROPERTIES OF THE 
QUASI-HYPERBOLIC MODEL 

From now on we consider only the case A = 0, e = —1 
and only expanding models. The corresponding conclu- 
sions for collapsing models follow immediately. 

It was stated in Ref. [l| that the surfaces H2 of con- 
stant t and z in (|2.8[) in the quasi-hyperbolic case e = — 1 
consist of two disjoint sheets. This was a conclusion from 
the fact that with e = — 1 the equation £ = has a so- 
lution for (x, y) at every value of z, and every curve that 
goes into the set 8 = has infinite length. However, it 
will be shown in Sec. IVII that the two sheets are in fact 



two coverings of the same surface, also in the general non- 
symmetric case. Their spurious isolation is a property of 
the stereographic coordinates used in (|2.8p . 

Note that with e = — 1 (|2.8p shows that for the signa- 
ture to be the physical (H ), 



E(z) > 1/2, 



(3.1) 



is necessary, with E = 1/2 being possible at isolated 
values of z, but not on open subsets. We shall also assume 



M(z) > 



(3.2) 



for all z, since with M < (|2.9p would imply i?, tt > 0, 
i.e. decelerated collapse or accelerated expansion, which 
means gravitational repulsion. 

In consequence of (|3.ip . only one class of solutions of 
is possible in the quasi-hyperbolic case: 



R 

t-t B 



M 
2E 



(cosh 77 — 1), 



M 



(2^)3/2 



(sinh?7 — 77). 



(3.3) 



The second of the above determines rj as a function of t, 
with z being an arbitrary parameter, and then the first 
equation determines R(t,z). 

Equations (pTTj) , and ([23)1 with A = imply that 
R,t 2 > at all z, i.e. there can be no location in the 
manifold at which _R, t = 0. In particular, there exists no 
location at which R — permanently. The function R 
attains the value only at t = is, i-e. at the Big Bang. 
We have, at all points where M > 0, 



lim R(t, z) 
t^t B 



0, 



lim R, t (t, z) 

t—¥tn 



00, 



(3.4) 



but R > at all values of z where t > ts- Thus, in 
the quasi-hyperbolic model there exists no analogue of 
the origin of the quasi-spherical model or of the center 
of symmetry of the spherically symmetric model. (This 
fact was demonstrated in Ref. [l| by a different method.) 

However, a location z = z m o at which M(z m o) = 
is not prohibited, even though the parameter 77 in (|3 .3[) 
becomes undetermined when M = / £. Writing the 
solution of (|3.3D as 



t-t B 



M 



y/4E 2 R 2 /M 2 + 4ER/M (3.5) 



(2£)3/2 

In (2ER/M + 1 + V / 4E 2 R 2 /M 2 + AER/M 



(where the log-term is the function inverse to cosh) we 
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sec that the limit of this as M — > is 3 



R 



lim (t — tn) = 



(3.6) 



and the same result follows from (12.91) with M = = A 



Note that (13. 6|) implies i?, t (t, z mo ) = ^/2£J (z mo ) - an 
expansion rate independent of time. This agrees with 
Newtonian intuition - expansion under the influence of 
zero mass should proceed with zero acceleration. At all 
other locations, where M(z) > 0, the expansion rate 
is greater than y / 2E(z), and tends to y/2E(z) only at 
R{t, z) — > oo. However, E(z) at z ^ z mo may be smaller 
than E(z mo ), so the expansion rate in the neighbourhood 
of the M = set may in fact be smaller than y/2E(z mo ). 

Conversely, at a location where R lt — constant (with 
A = 0), (123|) implies that M = (because i?, t 7^ in 
consequence of E > 1/2, M > and R > 0). 

The set where M = may or may not exist in a given 
quasi-hyperbolic Szekeres spacetime. It should be noted 
that, if it exists, it is a 3-dimensional hypersurface in 
spacetime; unlike the origin in the quasi-spherical mod- 
els. The latter is a 2-dimensional surface in spacetime 
and a single point in each space of constant t because in 
the quasi-spherical case M — implies R — via the 
regularity conditions. 



For an expanding model R, t > 0, and only past- 
trapped surfaces can possibly exist, for which k^;^ > 0. 
Apart from shell crossings the first factor in (|4.1[) is pos- 
itive everywhere. Hence, (|4.1[) implies 



R. 



y/2E-l 



+ e > 0. 



(4.2) 



For e = +1, and with R, t > that we now consider, this 
is fulfilled everywhere. For e = — 1 we get 



R,t >2E-1>0 



(4.3) 



(the last inequality from IpTTjl ). With A = 0, (|4~3]) is also 
guaranteed to hold everywhere, by (|2.9[) , since M > 
and R > 0. This means that every surface of constant t 
and z in an expanding model is past-trapped at all of its 
points. But then, every point of the Szekeres manifold 
lies within one such surface. This, in turn, means that 
every point of the Szekeres manifold is within a past- 
trapped region. Therefore, the whole quasi-hyperbolic 
expanding Szekeres manifold is within a past apparent 
horizon. 

The fact of being globally trapped is a serious limi- 
tation on the possible astrophysical applications of the 
quasi-hyperbolic model. 



IV. NO APPARENT HORIZONS 

In Ref. Q it was shown that a collapsing quasi- 
hyperbolic Szekeres manifold is all contained within the 
future apparent horizon, i.e. that it represents the inte- 
rior of a black hole. This is consistent with the fact that 
the corresponding vacuum solution (the hyperbolically 
symmetric counterpart of the Schwarzschild solution) has 
no event horizons and is globally nonstatic [l|. 

Here, we consider expanding models, and an adden- 
dum is needed to the result reported above. Consider a 
surface of constant t and z in (|2.8[) , and a family of null 
geodesies intersecting it orthogonally. As shown in Refs. 
Q and jl6j |. the expansion scalar for this family is 



& >jU — 2 



R, Z 
~R~ 



£ 



( R^t 



\V2E-1 



(4.1) 



where e = +1 for "outgoing" and e = — 1 for "ingoing" 
geodesies; 4 eq. (|4.ip was adapted to e = — 1. 



3 The result 113.61 1 shows that the argument used in deriving the 
regularity conditions at the center for the Lemaitre - Tolman 
model in Ref. [26| . and repeated in Sec. 18.4 of Ref. [13, was 
incorrect. The value of the parameter r\ need not be determined 
at the center. However, the resulting regularity conditions are 
correct because they can be derived in a different way. Once we 
know that R = 0, M = and R oc A/ 1 / 3 at the center, the 
behaviour of E at the center follows from eq. 112.911 . 

4 Since the surfaces of constant t and z are infinite, this labelling is 
purely conventional in this case, but the two families are distinct. 



V. INTERPRETATION OF THE 
COORDINATES OF (^8)) 



In order to understand the geometry of (|2.8[) . we begin 
with the hyperbolically symmetric subcase, P, 2 = Q, z = 
S, z = 0. It is most conveniently represented as 



ds 2 = dt 2 - 



Rjdz 2 



2E 



- - R 2 (dtf 2 + sinh 2 Mip 2 ) . (5.1) 



The two supposedly disjoint sheets of a constant- (t, z) 
surface, in the coordinates of (12.81). are 



Sheet 1 
Sheet 2 



x_P\ 2 , (y-Q s 1 
S 

x-P\ 2 (y-Q^ 1 



> 1, 



< 1. (5.2) 



S J \ S 

The transformation from Sheet 1 to (|5.1[) is 

{x,y) = (P,Q) + 5 , coth(z?/2)(cos^,sin(p), (5.3) 

while the transformation from Sheet 2 is 

(x,y) = (P,Q) + S tanh(i9/2) (cos ip, sin y), (5.4) 

This shows that the two sheets are in truth two coordi- 
nate coverings of the same surface. The direct coordinate 
transformation between the two sheets is the inversion 



r P n ^ S 2 (x>~P,y'-Q) 
[x - P,y - Q) = — 



{x 1 - P) 2 + (y' - Q) 



(5.5) 



■5 



The circle separating the two sheets, {x — P) 2 + (y — 
Q) 2 = S 2 , on which £ = 0, corresponds to $ — > ±00 
in the coordinates of (|5.1j) . The center of this circle, 
(x, y) = (P,Q), which is in Sheet 2, is mapped by (|5.4p 
to <d = 0. The infinity, (x - P) 2 + (y - Q) 2 -t 00, which 
is in Sheet 1, is mapped by (15 .3[) also to 1? = 0. These 
relations are illustrated in Fig. [TJ 

There is no reason to allow negative values of $ in (|5.1[) 
because, as both (|5 . 3[) and (|5.4[) show, the point of co- 
ordinates (— if) coincides with the point of coordinates 
(■#, (p + 7r), so the ranges <& € [0,+oo) and (p € [0,2tt) 
cover the whole (1?, y>) surface. 

Curves that go through $ = are seen from (|5.ip to 
have finite length. At $ — we have det (g a p) = 0, but 
the curvature scalars given in Appendix [A] do not depend 
on so i9 = is only a coordinate singularity. 

The geometry of the (x, y) surfaces in (|2.8I) is the same 
in the hyperbolically symmetric case and in the full non- 
symmetric case with £ , 2 7^ 0. In the (x, y) coordinates 
of (|2.8[) . z9 = corresponds to £ = —1 in Sheet 2, which 
is clearly not a singularity, and to £ — > 00 in Sheet 1. 
This seems to be a singularity in (|2.8[) , but the curvature 
scalars are not singular there, as shown in Appendix |B| 
Also the set £ = seems to be singular in (|2.8[) . but the 
same formulae in Appendix [Blshow that it is nonsingular. 
Hence, also in the general case there is no reason to treat 
these two sheets as disjoint - they are two coordinate 
coverings of the same surface. 



VI. GEOMETRY OF SUBSPACES IN THE 
HYPERBOLICALLY SYMMETRIC LIMIT 

A. Hypersurfaces of constant z 



This can be embedded in a flat 3-dimensional 
Minkowskian space with the metric 



by 



ds M = dT 2 - dX 2 - dY 2 



dt \ , 
— dR, 
dR J 

Y = i?sintf. 



(6.4) 




(6.5) 



The embedding (|6.5I) projects a point of coordinates 
(Rj'd) and points of coordinates (i?,$ + 27rn), where n 
is any integer, onto the same point of the Minkowskian 
space (|6.4[) . However, these points do not coincide in 
the spacetime (|5.1I) - the identification of (R, $) with 
(R, $ + 2ir) is not allowed because the transformation 
-d — s- i9 + 2n is not an isometry in (15.1[) . Thus, the surface 
with the metric (|6.3I) is covered by the mapping (|6.5[) an 
infinite number of times. This shows that a hyperboli- 
cally symmetric geometry is a rather exotic and compli- 
cated entity. We shall see this feature further on, while 
considering other surfaces. 

Using (12.91) with A = we can write 



dt 
dR 



R 



2ER + 2M 



(6.6) 



and then the integral in (|6 . 5|) can be calculated explicitly: 
FG M 



T 



2E E V / 2E{2E + 1) 
where D is a constant and 



^2EF + V2E + lGj +D, 
(6.7) 



A hypersurface z 
tensor: 



z x = constant has the curvature F =' ^(2E +1)R + 2M, G =' V2ER + 2M. (6.8) 



def 



-^0202 
^2323 



^0303/ sinh 2 tf 

^2. 



RR,t 



3 i? 2 .323 = R 2 sinh 2 d (1 - R, t 2 ) 



(6.1) 



where (x°,x 2 ,x 3 ) = (t,-d,ip). Consequently, it is fiat 
when R, t = ±1 and curved in every other case (also when 
R,t— constant ^ ±1)- However, (|2.9[) implies that with 
R,t = ±1 we have M — and E — 1/2 (recall: we con- 
sider only the case A = 0). Such a subset in spacetime (if 
it exists) is a special case of a neck - see the explanation 
to Fig. [TU] later in this section. 

The metric of a general hypersurface of constant z is 



ds z 2 = dt 2 - R 2 {t, zx) (di? 2 + sinh 2 tfd<p 2 



(6.2) 



To gain insight into its geometry, we first consider its sub- 
space given by ip — ipo = constant. The z\ is a constant 
parameter within R and will be omitted in the formulae 
below. The corresponding 2-dimensional metric is 



ds Zl ^ dt" 



-R 2 {t)&d 2 



dt 
dR 



dR 2 -R 2 dd 2 . (6.3) 



The constant D can be chosen so that T = at R = 0. 
Figure [2] shows the graph of the surface given by the para- 
metric equations (|6.5p as embedded in the 3-dimensional 
space with the metric (|6.4I) . It is not exactly a cone, the 
curves T(i?) do have nonzero curvature 



d 2 T 



M 



dR 2 FG 3 ' 



(6.9) 



but it is so small everywhere that it would not show up 
in a graph. Note that the vertex angle of this conical 
surface is everywhere larger than 7r/4, since dT ' /dR > 1 
from (|6.5p . 

Suppose that we made T unique by choosing D as in- 
dicated under (|6.8p . The vertex of the conical surface in 
Fig. [2]corresponds to the Big Bang. If we want the image 
in this figure to correspond to the history of the Universe 
from the Big Bang up to now, then the upper edge of 
the funnel should be at T(i? p ), where R p corresponds to 
the present moment. But this R p depends on the value 
of z = Z\. Consequently, the height of the funnel will be 
different at different values of z. 
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FIG. 1: Relations between the (i?, ip)- and (x,y) maps of a constant-(t, z) surface in 1 12, 81 and ( 15. Il l ■ The arrow marked by I corresponds 
to the transformation 115.41 1 that maps the set # = to (x, y) = (P, Q). Arrow II shows that both 1 15. 31 1 and 115.41 1 map t? -> oo to the circle 
£ = 0. Arrow III corresponds to 1 15. 31 1 that maps $ = to £ — > oo. 




FIG. 2: The surface of constant z — Zi and constant ip = tpo 
in a spacetime with the metric 115.11 1. The embedding is in a 
Minkowskian 3-space with the metric 1 16.411 . The vertex at R = 
lies at the Big Bang. The circles represent the surfaces of constant t 
and 2 in 15. H . This embedding is not a one-to-one representation, 
the surface in the figure is covered with that of 16.31 1 an infinite 
number of times - see explanation in the text. 



Now we go back to (|6.2[) and consider a surface of con- 
stant i? = r?o- Writing Co = sinhi9 we can write the 
2-metric as 



ds 



At 
dR 



dR 2 -d (C R) 2 -(C R) 2 dep 2 , 
(6.10) 



and then the embedding equations are 

T = 





R 



X 



2EC 2 R + 2MC 2 
CoRcosip, Y = CoRsinip. 



dR, 



3.11) 



Now there is no multiple covering because tp is a cyclic 
coordinate also in spacetime, and the surface given by 
(16. lip looks qualitatively similar to that in Fig. [2l except 
that the presence of Co introduces some flexibility. The 
second line of (|6. 1 1|) shows that the explicit expression 
for T is (|6.7[) multiplied by Co, with (M,E) replaced by 
Cq 2 {M,E). The radius of a circle of constant R is now 
(CqR). The value of Co is any in (— oo,+oo). When 
Co — > 0, the surface degenerates to the straight line X — 
Y = 0. In order that CoT in (16.71) allows a well-defined 
limit Cq —> 0, the constant D must have the form 



D 



M In C 



+ 



C EJ2E (2EC 2 



Co' 



(6.12) 



where D\ is another constant, 
so that C T = at R = 0. 
From (|6.11l) we find 



Again, it may be chosen 



lim 

Co— >cx 



dT 



d(C R) 



1. 



(6.13) 



so in the limit Co — > oo the surface (|6.1ip becomes ex- 
actly a cone with the vertex angle 7r/4. However, with 
Co — > oo the whole cone recedes to infinity, as can be 
seen from (I6.7[) and (|6.11l) : the vertex of the cone, which 
is at R — 0, has the property limc7 ^oo {CqT)\ r=q = oo, 
even with the value of D corrected as in (|6.12[) . 

Note that the image in Fig. [5] will not change qualita- 
tively when we go over from the hyperbolically symmct- 
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ric subcase (|5.ip to the general (nonsymmetric) quasi- 
hyperbolic case (|2.7I) - (|2.8I) . Each hypersurface of con- 
stant z in it is axially symmetric, and its metric can be 
transformed to the form (|6.2[) . Moreover, any surface of 
constant z and y can have its metric transformed to the 
form (|6.3|) . The only change with respect to Fig. [2] is 
that the cone-like surfaces, while still being axially sym- 
metric, can have their vertex angles different at different 
values of z. 

For completeness, we now consider the special flat hy- 
persurface with R,t= 1, M = and E — 1/2 men- 
tioned below (|6.1|) . It can be all transformed to the 3- 
dimensional Minkowski form. The transformation to the 
Minkowski coordinates (r, X, Y) is 

t = i?coshz9, X = R sinh $ cos <p, 
Y = R sinh i9 sin </>. (6.14) 

The surfaces of constant i? are given by the equation 

r 2 - X 2 - Y 2 = R 2 . (6.15) 

These are two-sheeted hyperboloids when R > and a 
cone when R — 0. They intersect the r-axis horizontally, 
and all tend asymptotically to the cone R = as X 2 + 
Y 2 -> oo (see Fig. [3}. 




FIG. 3: An axial cross-section through the family of hyperboloids 
given by 1 16. 151 . 

The surface (p = of (|6.14|) is depicted in Fig. Note, 
however, that Figs. [3] and 0] are graphs of a Lorentzian 
space mapped into an Euclidean space, so geometrical 
relations of (|6.2j) are not faithfully represented. 

B. The R(t,z) curves 

Where M > 0, we have R, tt < from ([23]) with A = 0. 
Consequently, R as a function of t must be concave. The 
slopes of the curves R(t, z) at various z depend on E(z), 




FIG. 4: The subspace {z = constant,^ = 0} of the spacetime 
I l5.lt with R = t. The hyperbolae are curves of different constant 
values of R, the straight lines emanating from the origin are lines 
of constant i9. 



and their initial points at t = ts arc determined by is(z), 
so both can vary arbitrarily when we proceed from one 
value of z to another. Fig. [5] shows a 3-d graph of an 
example of a family of R(t,z) curves corresponding to 
different values of z. 




t = t B (z) 

FIG. 5: An exemplary collection of the R(t, z) curves for vari- 
ous fixed values of z. The bang time curve t = tg(z) must be a 
decreasing function of z to avoid shell crossings [l| . The rightmost 
line is straight, corresponding to M = 0. The other t(R) functions 
in this figure are given by 113.511 with tg(z) = 2 — 0.05z 2 . They have 
M increasing with z (so \R,tt | is increasing as a function of z) and 
E increasing with z (so R,t is increasing). The values of (M,E) 
increase from (0.0,0.6) in increments of (0.1,0.1). Note that all 
curves except the M = one hit the t = ts set with R,t — > oo. 



C. Hypersurfaces of constant t 

Formulae for the curvature of the spaces of constant t 
in (|2.7[) - (|2.8[) are (from Ref. [l[ , in notation adapted to 
that used here): 

3 n 3d 

-Kl212 — ^1313 
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R(R, z -R£, z /£)(E, z -2E£, z /£) 



(2E - l)£ 2 



''R 



2323 



2ER 2 



(6.16) 



where the coordinates are labelled as (x 1 , x 2 , x 3 ) = 
(z, x, y). Equations (|6.16|) show that a space of constant t 
becomes flat when E = 0, but then it has the Lorentzian 

signature (H ). Consequently, with the Euclidean 

signature, these spaces can never be flat. 

Now let us consider the surfaces H2 of constant t = to 
and ip = tpo in (|5.ip . For the beginning we will assume 
that R, z > for all values of z in the region under inves- 
tigation. Then we can write the metric of H2 as follows: 



ds 2 2 = 



[dR(t , 



2E- 1 



R 2 M 2 



(6.17) 



When E = 1, this is the metric of the Euclidean plane 
in polar coordinates (R, $). With some other constant 
values of E, this will be the metric of a cone (see be- 
low). With other functional forms of E, it is the metric 
of a rotationally symmetric curved surface on which d 
is the polar angular coordinate. We encounter here the 
same phenomenon that was described in connection with 
Fig. [5] in each case, a point of coordinates (R, 1?) and 
points of coordinates (i?, + 27m), where n is any inte- 
ger, are projected onto the same point of the plane, cone 
or curved surface, respectively. However, as before, these 
points do not coincide in the spacetime (|5.1[) . Examples 
of embeddings of H2 in the Euclidean E 3 are illustrated 
in Fig. |Hl 

Note that with E = constant 7^ 1 other interesting 
geometries come up. If we interpret $ as a polar coor- 
dinate, then the ratio of a circumference of a circle R = 
constant to its radius is 2ir\ / 2E — 1, which means that 
the surfaces (|6.17p are ordinary cones when 1/2 < E < 1 , 
and cone-like surfaces that cannot be embedded in a Eu- 
clidean space when E > 1. With E being a function 
of z, the cones and/or cone-like surfaces are tangent to 
the (z, $) surfaces at the appropriate values of z. Thus, 
with E > 1, the (2,1?) surfaces cannot be embedded in 
a Euclidean space. Whether this surface looks like the 
smooth surface of revolution in the lower panel of Fig. [6j 
or like a cone, depends on the behaviour of E{z) in the 
neighbourhood of the axis R = 0. But attention: if the 
value of to under consideration is such that to > ts{z) 
for all z, then the set R = is not contained in the space 
t = to- We will come back to this below. We can discuss 
the embedding when we write the metric (16.171) as 



ds 2 2 



2(1 -E) 
2E- 1 



d-R 2 + di? 2 -t-i? 2 di? 2 . 



3.18) 



Now it is seen that with 1/2 < E < lwe can embed this 
surface in the Euclidean space with the metric d«3 2 = 
dX 2 + dY 2 + dZ 2 by 



X = Rcosd, y = i?sintf, Z = ± 



2(1 — E) 



2E- 1 



di?, 





'X 



FIG. 6: Upper graph: A surface Hi of constant t and ip in the 
metric 115.111 in the case when E = 1. It is locally isometric to the 
Euclidean plane, but points having the coordinates (R, 1? + 2irn) do 
not coincide with the point of coordinates (R, -d), so the projection 
of H2 covers the Euclidean plane multiply. The multiple covering 
is depicted schematically. Lower graph: When E is not constant, 
the embedding of a surface of constant t and ip in the Euclidean 
space is locally isometric to a curved surface of revolution, with a 
similar multiple covering, also shown schematically. The surface in 
the figure is the paraboloid Z = R? that results when E(z(R)) = 
(2R 2 + 1)/(4_R 2 + 1), where z(R) is the inverse function to R(t , z). 



while with E > 1 we can embed it in the Minkowskian 
space with the metric ds3 2 = — dT 2 + dA" 2 + dy 2 by 



X = Rcos'& 1 r = i?sintf, T = ± 



2(E-1) 



di?, 



(6.19) 



2E- 1 

(6.20) 

For later reference let us note, from (I6.19[) , that 
dZ/dR ->■ when E -> 1 and \dZ/dR\ -> 00 when 
E — > 1/2. This observation will be useful in drawing 
graphs and interpreting them. 

The surfaces on which t and $ are constant look similar 
to the surfaces described above, with two differences: 

1. The coordinate ip changes from to 2ir also in the 
spacetime (I5.1|) . so there is no multiple covering of the 
surfaces in the Euclidean space. 

2. The circumference to radius ratio is this time 
27rsinh^V2i? — 1, so the transition from cones to cone- 
like surfaces occurs at sinhi? = \ j\/2E — 1. 

Now let us recall what was said in the paragraph con- 
taining (I3.4[) : R(t,z) becomes zero only at t = ts- At 
any t > ts, R > for all z, even at M = as (|3.6[) 
shows. Thus, the surfaces in Fig. [5] can extend down to 
the axis R — only if, at the given instant t = t\, the 
function is(z) attains the value t\ at some z = z\\ then 
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t = ts at z = zi, so R(ti,Z\) = 0. This is illustrated in 
Fig. [7] If ti > ts(z) at all z, then R(t2,z) is nowhere 
zero. Let i?o > be the smallest lower bound of R(t2, z); 
then R(t2,z) > Rq > at all z, and the surface shown 
in Fig. [5] has a hole of radius Rq around the axis. Since 
R,t > 0, Rq is an increasing function of t, and the radius 
of the hole increases with t. This is illustrated in Fig. [8] 





k 1 t-e 




~~^tB(z) 

\ t = t1 




\ z 




z1 \ ' 



FIG. 7: The hypersurface t = t\ has a nonempty intersection 
with the Big Bang set t = tn(z). The function ll(tj , z) attains the 
value &l z z\ . and the surface from Fig. [6] extends down to 
the axis R = 0. At t = t2 we have t > tg at all values of z, so 
i?(t2,z) is nowhere zero and has a smallest lower bound Rq > 0. 
This means that the corresponding surface from Fig. [6] will have a 
hole of radius Ro around the axis. Since R,t > 0, the radius of the 
hole increases with time, as shown in Fig. [8] 




FIG. 8: The surface of constant t and <p of l|5.1jl . from the bottom 
graph in Fig. \E\ depicted at two instants > tg (bottom graph) 
and t3 > t2 (top graph). The multiple covering of the paraboloid 
is no longer taken into account. The hole around the axis expands 
along with the whole surface. 

So far, we have considered R as an independent vari- 
able within the space t = to- Since it is a function of z, 



the parameter along the radial direction in Figs. [5] and 
[8] is in fact the coordinate z. Now let us recall that R is 
also a function of t and that at every z there exists such 
a t (t = £b), at which R = 0. Thus, as we consider the 
spaces t = to at consecutive values of to, the surfaces de- 
picted in those figures get gradually "unglued" from the 
Big Bang set (which is represented by the axis of sym- 
metry R = 0), and expand sideways. At the moment, 
at which to begins to obey to > ts for all z, the surface 
becomes completely detached from the axis and contin- 
ues to expand sideways. This is when the hole mentioned 
above first appears. 

Let us also note the double sign in the definition of 
Z, (|6.19p . which was not taken into account in Figs. 
H] and |S] It means that each of those surfaces has its 
mirror-image attached at the bottom. In summary, the 
evolution of those surfaces progresses as shown in Fig. 
IHl The functions used for this picture are M = 10|z| 3 , 
E = 0.6 + 0. 5e"l z l, t B = -10 3 |z| + 100. The time instants 
are (t 1 ,...,t 6 ) = (1, 50, 100, 300, 500, 700). 




Y 



FIG. 9: Evolution of the surface from the lower panel of Fig. [6] 
The figure shows the axial cross-section of the surface at several 
time instants, ti < ••• < te- The Big Bang goes off along the 
Z axis, beginning at the top and at the bottom, and progressing 
toward the middle. The instant corresponds to the last moment 
when the surface has no hole. Multiple covering not shown. 

The nondifferentiable cusp at the plane of symmetry 
is a consequence of the assumption R, z > 0: to avoid a 
singularity in the metric (|5.1[) , E > 1/2 must hold every- 
where, and then (|6.19j) implies |dZ/di?| < oo everywhere. 
This means that the upper half of the surface cannot go 
over smoothly into the lower half. 

Let us now consider the case when R, z = at some z = 
z n . To prevent a shell crossing at z n , E(z n ) — 1/2 must 
also hold, so that lim 2 ^ 2?i (R, z /\/2E — l) is finite. This 
implies that R, z | =0 for all t (i.e. that the extremum 
of R is comoving), and then M, z — E, z = at z = z n 
from (|2.9[) . This is an analogue of a neck - an entity well 
known from studies of the Lemaitre - Tolman model 



But, as remarked under (|6.20l) . we have dZ/dR — > ±00 
where E — > 1/2. The evolution then looks like in Fig. 
ITUl The functions used for drawing it are M = 10 2 |z| 3 , 
E = 0.5 + |z| 3/2 , t B = -10 3 |z| 2 + 100, and the time 
instants are {h, . . . , i 6 ) = (1, 50, 100, 200, 300, 400). 



Z 

t, t 2 t 3 u t 5 t 6 




t 6 



Y 

FIG. 10: The analogue of Fig. [9]for the situation when R, z = at 
some z = z n (in the middle horizontal plane) . Then the upper half 
of each constant-(t, ip) surface goes over smoothly into the lower 
half. 

The minimum of R with respect to z need not exist in 
any space of constant t. (By minimum we mean not only 
a differentiable minimum similar to the one in Fig. I1Q[ 
but also a cusp at the minimal value like the one in Fig. 

It will not exist when the function ts{z) has no upper 
bound, i.e. when the Big Bang keeps going off for ever, 
moving to ever new locations. In that case, the surfaces 
shown in the upper half of Fig. |H]will never get detached 
from the Big Bang set, only the vertex of each conical 
surface will keep proceeding along the axis. In the special 
case ts = constant, the whole surface of constant t and ip 
gets "unglued" from the axis R — at the same instant. 
The image would look similar to Fig. [TU1 but there would 
be no conical surfaces with the vertices progressing along 
R = 0. The generators of the surface in the picture are in 
general not vertical. They become vertical when R, z = 
everywhere, i.e. when R = R(t), which can happen only 
in the /3, 2 = family of Szekeres solutions that we do not 
discuss here. 



VII. SPACES OF CONSTANT t IN THE 
GENERAL QUASI-HYPERBOLIC CASE 

The main difference between the hyperbolically sym- 
metric case, where £, z = 0, and the full quasi- hyperbolic 
case, where £, z ^ 0, is seen in (|2.8[) . Consider two sur- 
faces Si and 5*2 such that t = to = constant on both, 
z = z\ on Si and z = Z2 on S^- When £ , z = 0, the 
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FIG. 11: Left: The view from above of the surface from the 
lower panel in Fig. \E\ The circles are images of the curves of 
constant R (and thus of constant z). Right: An example of a 
corresponding image in the general quasi-hyperbolic case. Now the 
geodesic distance between the circles depends on the position along 
the circle. 



geodesic distance between S\ and S2 along a curve of 
constant (x 7 y) is the same for any (x, y). When £, z ^ 0, 
this distance depends on (x,y) and varies as the func- 
tions P(z), Q{z) and S(z) dictate. Figure fTTl left panel, 
shows an exemplary family of constant-i? curves in a sin- 
gle surface of constant t and ip with £, z = 0. This is a 
surface of constant t and (y/x) in the coordinates of (|2.8[) 
- a contour map of the surface from the lower panel in 
Fig. |U With a general £{x,y,z) the geodesic distance 
between the constant-i? curves will depend on the posi- 
tion along each curve, and the whole family would look 
like in the lower panel of Fig. [TT] 5 



VIII. INTERPRETATION OF THE MASS 
FUNCTIONS M(z) AND M(z) IN THE 
QUASI-SPHERICAL CASE 

In the quasi-spherical case, the function M(z) of (|2.9|) . 
by analogy with the Newtonian and the Lemaitre - Tol- 
man cases, is understood as the active gravitational mass 
inside the sphere of coordinate radius z. In fact, it is puz- 
zling why it depends only on z when the mass-density 
(|2.10p so prominently depends also on t, x and y. Some- 
what miraculously, as shown below, the denominator in 
(|2.10[) is canceled by the \J—g\\ term inside the inte- 
gral J P\f\gz\d?,x that determines the mass in a sphere. 6 



5 Fig. 1111 is in fact deceiving. The curves shown there as circles 
are images of infinite curves, as explained under 116.511 . and each 
image is covered an infinite number of times. 

6 93 is the determinant of the metric of the 3-space t = constant 
in |gl}. 
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The term containing £ , z in the numerator gives a zero 
contribution to the integral. This is consistent with the 
fact, known from electrodynamics, that the total charge 
of a dipole is zero (see Refs. [U, [l?], HqJ for the split- 
ting of (|2.10l) into the monopole and the dipole part in 
the quasi-spherical case). It is also consistent with the 
result of Bonnor [§, |i| that the Szekeres solution can be 
matched to the Schwarzschild solution. 

The considerations of this and the next three sections 
arc intended to prepare the ground for an analogous in- 
vestigation in the quasi-hyperbolic case further on. The 
questions we seek to answer are: Can M still be inter- 
preted as mass, and where does the mass M{z) reside 
when a surface of constant z has infinite surface area? 

Let us calculate, in the quasi-spherical case, the 
amount of rest mass within the sphere of coordinate ra- 
dius z at coordinate time t, assuming that z = zq is 
the center, where the sphere has zero geometrical radius 
(see Ref. [la]). This amount equals M = § v p^J\g?~\d3X, 
where V is the volume of the sphere and p is the mass 
density given by (|2.10[) . Substituting for p and g% we get 



M 



+ 00 



— / dx I dy I du 

^ n J — OO J — OO J Zq 

M, u («) 3M£ ,„ 



VI + 2E£ 2 VT+2ES 



(8.1) 



where u is the running value of z under the integral. Note 
that £ is the only quantity that depends on x and y, it is 
an explicitly given function, and so the integration over 
x and y can be carried out: 



IX. THE MASS IN THE SPHERICALLY 
SYMMETRIC CASE 

In nearly all the papers concerning the LT model and 
the quasi-spherical Szekeres model it was assumed that 
each space of constant t has its center of symmetry (in the 
LT case) or origin (in the quasi-spherical Szekeres case), 
where M = and R = at all times. We will assume the 
same here, but this is an assumption. It is possible that 
the center of symmetry is not within the spacetime in the 
LT case, and the corresponding quasi-spherical Szekeres 
generalization will then have no "origin" . 8 

For the beginning we will consider the spherically sym- 
metric (Lemaitre-Tolman) subcase, in which P lZ = Q, z — 
S, z = 0, so P = Q = can be achieved by coordinate 
transformations. Then £, z = in (|8.2[) . Now suppose 
that we calculate the integral in the first of (18.21) over 
a circular patch C of the sphere (circular in order that 
no (x,y) dependence appears from the boundary shape). 
The boundary of C is an intersection of the sphere with 
a cone whose vertex is at the center of the sphere. Let 
the vertex angle 9 of the cone be 7r/n. This translates 
to the radius of C in the original (x, y)-coordinates being 

i*o = S'tan(7r/2n) *= Sf3. Then we have in place of the 
first of (El 



= 2tt [1 — cos(7r/n)] 



= 47T- 



/3 2 



1 + /3 2 



(9.1) 



-OO r + OO £ 

dx / dy-^f = 0. 

OO J — OO J — OO J — OO ^ 

(8.2) 

(The first of these just confirms that this is the surface 
area of a unit sphere.) Using this in (|8.1|) we get 



M = 



M„ 



V1 + 2E 



(u)du, 



(8.3) 



which is the same relation as in the LT model, 7 and shows 
that l/y/l + 2E is the relativistic energy defect/excess 
function (when 2E < and 2E > respectively). 

In the quasi-spherical case we are able to calculate the 
integral with respect to x and y over the whole (x,y) 
surface because its surface area is finite. Such a calcula- 
tion cannot be repeated for the e < cases because both 
integrals analogous to (|8.2[) are infinite. Let us then con- 
sider what happens with M and A4 when we calculate 
the integrals in (|8.2I) over a part of the sphere. 



This tends to 47r when n — > 1 (uo — > 00). Note that 
the final result does not depend on S - this happened 
because we have chosen the coordinate radius in each 
circle, \/ x 1 + y 2 = u , to be a fixed multiple of S. 

In the spherically symmetric case now considered, we 
choose the same cone to define the circles of integration 
in (|9.1|) in all surfaces of constant z. Instead of (18.31) we 
get for the amount of rest mass within the cone, A4c : 



M c = 



du- 



1 



M„ 



■ /3 2 J zo VT+2E 



(u)du. 



■2) 



The term of J8JJ) that contained £ lU disappeared here in 
consequence of the assumed spherical symmetry, but it 
will not disappear when we go over to the general case, 
and its contribution will have to be interpreted. 



For the quasi-spherical case also other integral relations are sim- 
ilar to the ones in the L-T model, which is a consequence of the 
fact that the dipole contribution vanishes after averaging over 
x-y surface [27]]. 



The authors are aware of just one paper, in which LT models 
without a center of symmetry were considered. These are the 
"in one ear and out the other" and the "string of beads" models 
of Hellaby [28l |. described also in Ref. [l7| . In both of them, 
M ^ throughout the space. 
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X. SYMMETRY TRANSFORMATION S OF A 
SPHERE IN THE COORDINATES OF <f2TT2]> 

In the spherically symmetric case rotations around a 
point are symmetries of the space. In this case, if we 
rotate the whole cone around its vertex to any other po- 
sition, eq. (|9.2|) will not change. Each z — const circle 
will then be rotated by the same angles to its new po- 
sition, and the result of such a rotation will be a cone 
isometric to the original one. 

However, in the general, nonsymmetric case the 
spheres z = constant are not concentric. Suppose we 
build, in the general surface composed of circles, 

each circle taken from a different sphere. If we rotate 
each sphere by the same angles, whatever surface existed 
initially, will be deformed into a shape non-isometric to 
the original one. We want to calculate the effect of such a 
transformation, and for this purpose we need the formu- 
lae for the 0(3) rotations in the (x, y)-coordinates. We 
calculate them now. 

The generators of spherical symmetry, in the ordinary 
spherical coordinates, are (l7| : 



d d d 

J\ = J2 =sm<p— +cosy?cot#— , 

dip ov o<p 



9 a 9 
J3 = COS W— — Sin W COt I? Tr- 



(10.1) 



We transform these to the (x, y)-coordinates of (|2.12l) . 
for the beginning with P = Q = 0, S = 1, by 

x = cot(i?/2) cos(^), y = cot(i?/2) sin(^) (10.2) 

and obtain 

d d 
■h = x— - y—, 
ay ox 

9 1 9 9\ d 
J-2 = 2xy— + (l-x 2 +y 2 )—, 
ox 'ay 

d d 

J 3 = (l + x 2 -y 2 )- + 2xy-. (10.3) 

The transformations generated by J\ are rotations in the 
[x, y)-plane. To find the transformations generated by 
J3 we have to solve (see Ref. [13] for explanations): 

% = l + X ' 2 - y ' 2 > % = W> (1 °- 4) 

where A is the parameter of the group generated by J3. 
The general solution of this is 



y 

x' 



l/U; 



Ui d = C + \/C 2 - lcos(2A + D), 
VC 2 - 1 sin(2A + D)] /U x , (10.5) 



where C and D are arbitrary constants of integration. 
These have to obey the initial conditions (a?', 2/)L=0 = 
(x,y). After solving for C and D this leads to 



c' = [2xcos(2A) + (1 



y I sin 



(2A)] /U 3 



U 3 d = 1 + x 2 + y 2 + (l - x 2 - y 2 ) cos(2A) 
-2irsin(2A), 

V - 2y/U 3 . (10.6) 

It is instructive to calculate the effect of the transfor- 
mation (| 10 . 6[) in the (fl,<p) coordinates of (|10.2[) . Let 
us then take a point of coordinates (x,y) = (a;, 0), i.e. 
(■d,tp) = (i?,Q), and let us apply (TiTTo) to it. After a 
little trigonometry we get 



tan(i7'/2) 



1 — cos $ cos(2A) — sini?sin(2A) 



sint?cos(2A) 
tan(i?/2 - A) = 



cos"i9sin(2A) 

= i?-2A, (10.7) 



i.e. (|10.6|) is equivalent to rotating the sphere around the 
(i?,<p) = (7r/2,7r/2) axis by the angle (-2A). 
It can now be verified that the quantity: 



def 1 



V 



■ x' 2 + y' 2 



2y 



2y> 



(10. 



is an invariant of the transformation (|10.6j) . The set I = 
2C is the circle x 2 + (y - C) 2 = C 2 - 1. 

An arbitrary circle of radius A and center at x = y = 0, 



y 



A 2 , is transformed by (jl0.6l) into the circle 



'1 + A 2 ) sin(2A) 



1 2 



(l + ^l 2 )cos(2A) + l-A 2 

4^2 

[(l + A 2 )cos(2A) + l"^l 2 ] 2 ' 



(10.9) 



The coordinate radius of the (x',y r ) circle is different 
from the original radius A except when A = 0, which is 
the identity transformation. But the coordinate radius 
is not an invariantly defined quantity. An invariant mea- 
sure of the circle, its surface area, does not change under 
the transformation (|10.6|) , and neither does the invariant 
distance between any two points, as we show below. 
The Jacobian of the transformation (|10.6[) is 



d(x',y') 
d(x,y) 



U, 



(10.10) 



Together with (|10.8[) and (|10.6I) this shows that the sur- 
face element under the integral (19.21) . Adxdy/E 2 , does not 
change in form after the transformation. (This must be 
so, since (|10.6I) is just a change of variables that does not 
change the value of the integral.) Thus, (|10.6I) preserves 
the area of any circle, as is appropriate for a symmetry. 

The invariant distance between the points {x, y) 
ft). 0) and (x, y) = (A, 0) (i.e. the invariant radius of 
the original circle referred to in (|10.9[1 ) is, from (12.71) 
flU} with £ = +1, P = Q = 0, S = 1: 



dx 



2 arctan A. 



(10.11) 
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The image under (|10.6|) of any point (x, 0) is (xi(x),0), 
where, from (|10.6|) : 



xi(x) = 



2a;cos(2A) + (l - x 2 ) sin(2A) 



1 + x 2 + (1 - x 2 ) cos(2A) - 2x sin(2A) ' 

(10.12) 

Thus, the image of (0,0) is (xq,0), where 



x 



sin(2A) 
l + cos(2A)' 



(10.13) 



The invariant distance between the images of (0, 0) and 
of (A, 0) is then 



Xl(A) d Xl 

xi (0) 1+Zl 2 



2 arctan(xi)| 2 .J 



xi{A) 



2 arctan A, 



(10.14) 

by employing the identity arctan a — arctan /3 = 
arctan [(a — (3)/(l + a/3)]. This certifies that the invari- 
ant distance between the center of a circle and a point on 
the circle is the same as the invariant distance between 
their images (but the image of the center is no longer the 
center of the image-circle, compare (|10.9|) and (|10.13[) ). 

The transformations generated by J2 result from those 
for J3 by interchanging x' with y' and x with y; then 
all the conclusions about invariant properties follow also 



for these transformations, and, in consequence, for any 
composition of (|10.6[) with them. 

XI. THE MASS IN THE GENERAL 
QUASI-SPHERICAL CASE 

Now let us consider the general case, and integrals 
analogous to (|8.2[) . where the (x, y) integration extends 
only over a circular subset of each sphere, the radius of 
each circle being a fixed multiple of S. In the general 
case, each sphere has a geometrically preferred center at 
(x, y) — (P(z),Q(z)), and, for the beginning, we choose 
the center of the disc of integration C at that point. As 
before, the radius of each circle will be a fixed multiple 

def 



of S: \Jx 2 + y 1 



Uq 



S(3. This means, this time 



the volume of integration will not be a simple cone, but 
a 'wiggly cone' - the circles in the different z = const 
surfaces will have their centers not on a straight line or- 
thogonal to their planes, but on the curve given by the 
parametric equations x = P(z), y = Q(z) that is not 
orthogonal to the planes of the circles. The result (|9 . 1 [) 
still holds within each z = const surface, but the ana- 
logue of the second integral in (|8.ip . calculated over the 
interior of the 'wiggly cone' here, will no longer be zero. 
Instead, introducing in each z = const surface the polar 
coordinates x — P + u cos ip, y = Q + u sin ip, we get 



f £, s ( 2lx p (S, z /2)(l~u 2 /S 2 )- j,(P, z u cos p + Q, z u sin^) 
J duv/— / clip I ndii - 



8nSS, 



Jo 

)2 



u ° uS, z (S 2 -u 2 ^ 







(S 2 



-du = 4irSS, z 



(S 3 /8)(l + u 2 /S 2 ) 3 



Uo 



(S 2 + u 2 Y 



= 47T- 



S (I + /32) 



2 ■ 



(11.1) 



In agreement with (|8.2|) this goes to zero when uq — > 00. Consequently, from (|8.ip . (|9.2[) and (jll.l|) . the total mass 
within the wiggly cone is 



M = 



8 2 



1 + P 2 J ZQ 



M, u B 2 
(u)au — 6- 



z MS,- 



(l + /3 2 r Jzo SVT+2E 



[u)du. 



(11.2) 



r 



It contains a contribution from S, z that is decreasing with 
increasing /3, i.e. the greater volume we take, the less 
significant the contribution from S, z gets. It will vanish 
when the integrals extend over the whole infinite range 
of x and y (in the limit /3 — > 00). This can be interpreted 
so that in a wiggly cone the dipole components of mass 
distribution do contribute to M. - but less and less as the 
volume of the cone increases. Thus, with such choice of 
the integration volume M does not have an immediate 
interpretation - but it becomes proportional to the mass 
within the cone in the spherically symmetric limit. 

As an example, consider the axially symmetric family 
of spheres whose axial cross-section is shown in Fig. 
The circles are given by the equation 

(x-\/b 2 + u 2 ^ 2 + y 2 = u 2 , (11.3) 



where & is a constant that determines the center of the 
limiting circle of zero radius, while u is the radius of the 
circles. (The same family of spheres was used in Ref. [l[ 
to construct Szekeres coordinates for a flat space.) Figure 
[TBI left graph, shows the initial wiggly cone constructed 
for these spheres - the one referred to in (|11.1[) and (jll.2l) . 



We derived the transformation (|10.6p in the coordi- 
nates in which the constants (P, Q, S) were set to (0, 0, 1) 
by coordinate transformations. With general values of 
(P, Q, S), the result would be 
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x' - P 
S 



S 

y-Q 



■cos(2A) +11- 



(x - P) 2 + (y - Qf 
S 2 



sin(2A) 



/U 4 



TT dcf (x-P) 2 + (y-Q) 2 
1/4 = 1 + ^ 



1 - 



(.t - p) 2 + ( y - g) 5 



S 2 



cos(2A) - 2^-— -sin(2A). (11.4) 



Now let z = Z\ correspond to the base of the wig- 
gly cone, where the values of the arbitrary functions are 
Pi = P(zi), Qi = Q(zi) and Si = S(zi). Apply the 
transformation (|11.4j) with (P,Q,S) = (Pi , Q x , Si) to 
each sphere intersecting the wiggly cone. In the base 
z = Zi this will be a symmetry, in other spheres this will 
not be a symmetry. One should in principle calculate 
the effect of this transformation in other spheres on the 
integrands in (|9.2j) and (jll.lj) to see what happens. But 
then, (111 .4)) is merely a change of variables under the in- 
tegral that does not change the value of the calculated 
integral. Thus, eq. (jll.2p applies independently of where 
we choose the base of the wiggly cone, and its position 
that we chose initially (each circle had its center in the 
geometrically distinguished center of the (x, y)-surface) 
is only necessary to fix the relation between circles cor- 
responding to different values of z. 

The result of the transformation (jll.4l) applied to the 
wiggly cone of the left graph in Fig. [13] is shown in the 
same figure in the right graph. 
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FIG. 12: An axial cross-section through the family of spheres 
given by 111.31 . 



XII. THE SYMMETRY TRANSFORMATIONS 
IN THE HYPERBOLIC ALLY SYMMETRIC CASE 

Before proceeding to the case of interest, let us con- 
sider the hyperbolically symmetric subcase, in which 
Piz = Qiz = S, z = 0, and so P = Q = by a transfor- 
mation of x and y. The symmetry group of the resulting 
metric is a subgroup of that for the corresponding vac- 
uum solution [lj: 



ds 2 = - 1 



2™ \ , „o 
~RJ 



+ 2m/R 
- R 2 (d-d 2 + sinh 2 Mip 2 ) , 

The full set of Killing vectors for this metric is 

«i = do , fc 4 = d 3 , 

k2° — cos</5d~2 a — coth$sint£?53 a , 

k^ a = sin ipS2° + cothi? cos (fiS^ , 



dR 2 



(12.1) 



(12.2) 



where (x , x , x 2 , x 3 ) = (T, R,$, tp). To find the symme- 
tries explicitly, we have to transform (|1 2 .2|) to the coor- 
dinates of (j2"~g|) . 

We are interested in the transformations within an 
(x,y) surface, and those are generated by k2, &3 and k±. 
The transformation from the ($, <p) coordinates of (|12.1[) 
to Sheet 2 of the (x, ^-coordinates of (|2~%j) is (pT4"l) . The 
inverse formulae are 



cosh?9 = 
sinh $ = 



arctan(y/a;), 



1 - {x 2 + y 2 ) ' 

2\Jx 2 + y 2 
1 - (a; 2 + y 2 ) ' 



(12.3) 



We now transform the Killing fields (PUg) by (TT2~5)) . In 
the (x, y)-coordinates we get for the generators 



d d 
Ji = x— - y—, 
oy ox 



J2 = [1 



■h = 



ax 



- X ^dx- 2xV dy-> 



1 



(12.4) 
(12.5) 
(12.6) 



The Ji generates rotations in the (x, y) surface. To find 
the transformations generated by J2, we have to solve 
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FIG. 13: Left graph: a wiggly cone constructed for the family of spheres shown in Fig. 1121 The vertex angle for this cone is 7r/32. 
Right graph: the result of the transformation l|11.4|l applied to the cone from the left graph. The initial cone is shown in thin lines. The 
rest mass contained in the new wiggly cone is the same as it was in the initial cone. Dotted lines show the image of the initial cone that 
would result if each circle of intersection of the original cone with a sphere were rotated by the same angle around the center of its sphere. 



the set 

da/ 
dA 



1 + y — x , 



ay 

dA 



-2x'y' 



(12.7) 



with the initial condition that at A = we have (a/, y') 
(x,y). The general solution of this set is 



y' = 1/Wi, Wi = -C+ VC 2 + lcosh(2A + L>), 
x' = (l/VKiVC 2 + lsinh(2A + D), (12.8) 

where C and D are arbitrary constants to be determined 
from x'(0) = x, y'(0) = y. They are 

sinh-D = 2x/W 2 , cosh£> = (x 2 + y 2 + l) /W 2 , 

W 2 = ^(x 2 +y? - I) 2 +4y 2 , 

C = 



1 



2y 



(x 2 + y 2 -l). 



(12.9) 



So, finally 

x' = (1/Wa) [2a;cosh(2A) + (l + x 2 + y 2 ) sinh(2A)] , 

W 3 d = 1 - (x 2 + y 2 ) + (1 + x 2 + y 2 ) cosh(2A) 

+2a;sinh(2A), 
y' = 2y/W 3 . (12.10) 



The equations corresponding to (|12.7p for the genera- 
tor J3 result from (|12.7p simply by interchanging x' with 
y' . The corresponding initial condition then results by 
interchanging x with y. Thus, from (|12.10p we can read 
off the transformation generated by J3; it is 

x = 2x/Wi, 

W 4 d = 1 - (x 2 + y 2 ) + (1 + x 2 + y 2 ) cosh(2A) 



+2y sinh(2A), (12.11) 
y' = (I/W4) [2ycosh(2A) + (l + x 2 + y 2 ) sinh(2A)] . 



Note that the (x',y') given by (|12.10|) also obey the 
third equation in (|12.9I) . so the quantity 



7- def 1 / 2 

2y 



1) 



(12.12) 



is an invariant of the transformations (|12.10[) . The cor- 
responding invariant for (|12.11l) is 



T dot 1 / 2 

h = = — [x 
2x v 



y 



(12.13) 



These facts are helpful in calculations, and so is the fol- 
lowing identity that follows from ()12.12[) 



/2 



w 3 



(12.14) 



The fact that l\ is an invariant of (|12.10p means that 
the transformation (|12.10|) maps the set 1% = C = con- 
stant into itself for every C . This set is a circle of radius 
\JC 2 + 1 and center in the point (x, y) = (0, C). 

The inverse transformation to (I12.10P results from 
(|12.10l) by the substitution A — > (—A). This can be veri- 
fied using the above identities. 

These same identities can be used to verify that the 



transformation (|12.10|) maps the circle 



y 1 = A 2 into 



the following circle in the (a/, y')-coordinates 



(1 - A 2 ) sinh(2A) 



l + A 2 + (l- J 4 2 )cosh(2A) 

4^ 

[1 + A 2 + (1- A 2 )cosh(2A)] 2 ' 



(12.15) 
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The radius of this new circle equals the original radius, 
A, only in two cases: A = 0, which is an identity trans- 
formation, or A = 1. In both cases, also the center of 
the circle remains unchanged. The radius meant here is 
a coordinate radius that has no invariant meaning. The 
meaningful quantity is the geometric radius, which is the 
invariant distance between the center of the circle and a 
point on its circumference. It can be verified that the 
invariant distance between any pair of points is the same 
as the invariant distance between their images. 

Below we present some remarks about the transforma- 
tion (|12.10[) . The same statements apply to (|12.11j) . 

The Jacobian of the transformation (I12.10|) is 



d(x,y) 




(12.16) 



which, together with (112. 14|) . shows that the integrand in 
the integral J d 2 xy/£ 2 is form- invariant under this trans- 
formation. In this integral, (|12.10[) is an ordinary change 
of variables, and the area of integration in the (x',y')~ 
variables will be an image of the original area under the 
same transformation. This means that also the value of 
this integral is an invariant of (|12.10[) . Thus, if we choose 
the region of integration to be a circle around a point, it 
does not matter where the center of the circle is because 
we can freely move the circle around the (x, y) surface by 
symmetry transformations without changing its area. 



XIII. THE MASS FUNCTION IN THE 
QUASI-HYPERBOLIC CASE 

In the quasi-hyperbolic case, the {x, y) surfaces are in- 
finite, so they do not surround any finite volume. Thus, 



unlike in the quasi-spherical case, we should not expect 
the value of the mass function M(z) to correspond to 
a mass contained in a well-defined volume. We should 
rather observe the analogy to a solid cylinder of finite ra- 
dius and infinite length in Newton's theory, in which the 
mass density depends only on the radial coordinate. Its 
exterior gravitational potential is determined by a func- 
tion that has the dimension of mass, whose value is pro- 
portional to mass contained in a unit of length of the 
cylinder. 



We now proceed by the same plan as we did in the 
quasi-spherical case. We can freely move a circle of inte- 
gration around each (x, y)-surface. We first consider the 
hyperbolically symmetric case and we erect over a cho- 
sen circle a solid column in the z-direction that contains 
a certain amount of rest mass. Then we go over to the 
quasi-hyperbolic nonsymmetric space and erect a wiggly 
column that contains the same amount of rest mass. 



We will integrate over the interior of a circle in Sheet 2 
whose radius uq is, for the beginning, unknown. We only 
know that the radius must be smaller than 5, so that 
the integration region does not intersect the circle where 
£ = (since, we recall, £ — is the image of infinity, and 
the integral over a region that includes £ = would be 
infinite). Thus, instead of (|9.1[) and (jll.ip we have this 
time 



f £jz f * f U ° ~SuS 2 , , Or, „„„, „„ ,_\ 4ttSS, z Uq 2 

/ d 2 xy — = dip — —3 (u cos (pP, z +u sm ipQ, x +u S, Z /{2S) +SS, z /2) du = — — — j. 

iv c Jo Jo (u - &> z ) (wo - o ) 



The first integral in (| 13 . 1 1) will be independent of S when 
Mo is a fixed multiple of S: 



Then 



u = /3S, /3 < 1, 



d 2 xy — — 47r 



v 



£ 2 1-/3 2 ' 



(13.2) 



4n(S, z /S)f3 2 
= (1-P 2 f ' (13 - 3) 



Instead of (|11.2[) we now have: 



M = 



P 2 



M, u («) 



du- , 



du 



3MS, U 
Sy/2E - 1 ' 



(13.4) 



The meaning of the limits of integration has to be 
explained here. In the quasi-spherical case, and in the 
spherically symmetric LT subcase, one usually assumes 
that each space of constant t has its center of symmetry, 
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where M = = R. As explained at the beginning of Sec. 
IIX1 this is an additional assumption - the center of sym- 
metry need not belong to the spacetime. But the center 
of symmetry, or origin, is the natural reference point at 
which the mass function has zero value. In the quasi- 
hyperbolic case now considered, a similar role is played 
by the set M = 0, so we will assume that this set exists. 
Again, as mentioned earlier, this set is a 2-dimensional 
surface in each space of constant t, and not a single point. 

With this assumption made, zq in (|13.4[) will be the 
value of z at which M(zq) — 0. In addition, we assume 
that (M, u /V2E- 1) and [S, u /{S\/2E- 1)] are finite at 
z = zq and in a neighbourhood of zq, so that both inte- 
grals in ()13.4|) tend to zero as z — > z . 

These equations are very similar to the corresponding 
ones in the quasi-spherical case, (|11.1[) - (|1 1 .2|) . so one 
is tempted to interpret M , by analogy with that case, as 
a quantity proportional to the active gravitational mass 
contained within a solid (wiggly) tube with circular sec- 
tions, the radius of a circular section at z = Z\ being 
/3S(zi). The base of the tube is in the surface (t = to = 
const, z), its coordinate height is (z — zq), and its top 
is at (t, z) = (to,Zo). There is no problem with this in- 
terpretation in the hyperbolically symmetric case, where 
S, u — and the second integral in (|13.4|) vanishes. 

However, there is a significant difficulty when going 
over to the quasi-hyperbolic case. In the quasi-spherical 
case we were free to take the limit of the integral extend- 
ing over the whole sphere, which was /3 — > oo. Here, the 
integral over the whole hyperboloid would correspond to 
/3 1, and in this limit all the integrals (lT3~Tj) - ((TO)) 
become infinite. Worse still, the contribution from the 
dipole - the second integral in (|13.4I) - tends to infinity 
faster than the monopole component (the first integral), 
while in the spherical case increasing the area of integra- 
tion caused decreasing the influence of the dipole. 

We can do another operation on (113.41) that will shed 
some light on the meaning of M. The volume of the 
region containing the rest mass M. is 



V = / du d 2 xyy/\g 3 (t,u,x,y)\, (13.5) 
J z Ju 

where g% is the determinant of the metric of the 3- 
dimensional subspace t = constant of (|2.ip . thus 



(13.6) 



i , f i R 2 {R,u — R£ ,u /£ ) 



This has the same structure as the integral representing 
M. Since R and E do not depend on x and y, the inte- 
gration with respect to (x,y) can be carried out, and by 
(ITOi - (fT33)) we get 



V 



4tt^ 2 f z R 2 R 



du 



1 - P 2 J zo V2E—J 
4tt/3 2 f z R 3 S, U 



(1 - /3 2 ) 2 Jz SVZE^T 



du. (13.7) 



For the ratio M. /V we now calculate two consecutive lim- 
its: first /3 — > 1, to cover the whole of each z — constant 
hyperboloid, and then z — > z^, where z^ is the value 
of z at which R — > oo, to cover the whole t — constant 
space. After taking the first limit, we get 



lim 



M 
V 



3MS, U 
Sy/2E - 1 



■.du 



4tt 



=du 



sVTe- i 

(13.8) 

Since in general (apart from special forms of the functions 
involved) both the numerator and the denominator above 
become infinite when z — > z^, we apply the de l'Hopital 
rule and obtain 

M . 3M 
lim lim — = lim — — . (13.9) 

The l.h.s. of the above is the global average of rest mass 
M. per volume. The r.h.s. looks very much like the same 
type of global average for the active gravitational mass 
M, except that it is taken with respect to a flat 3-space. 

A very similar result follows when we take z — ¥ zq 
instead of z — > z^ in (|13.8[) . Then both the numerator 
and denominator tend to zero and we obtain 



lim lim 



M 
V 



3M 
lim -. 

z-^zq AnR fi 



(13.10) 



On the l.h.s here we have a global average of M. /V over 
the (x, y) surface taken at the value of z at which M = 0. 

Equation (|13.10p results also when the integrals in 
(]13.4[) - (113. 8[) are taken over the interval [z\ , z 2 ] , where 
zq < z\ < Z2 < Zoo, and then the limit z 2 — > z\ is taken. 

All the calculations in this section were done in Sheet 
2 of the (x, y) map. The corresponding results for Sheet 
1 are obtained by taking all integrals with respect to u 
over the interval [uo,oo) (with uq > S now) instead of 
[0,uo], substituting 1/uq for u in (|13.1[) . and 1//3 for (3 
(with the new /3 obeying (3 > 1) in P^|) . pXl]) and 
P^7|) . Equations pX5|) - (|13.10p do not change. 

The meaning of the limits on the r.h. sides of (|13.9[) 
and (| 13 . 1 0[) requires further investigation. Note that they 
arise from the dipole contributions to mass and volume. 



XIV. SUMMARY 

The aim of this paper was to clarify the geometrical 
structure of the quasi-hyperbolic Szekeres models [1, H| 
given by (|2.7[) - (|2.9p . and of the associated hyperbol- 
ically symmetric dust model given by (|5 . 1 [) . The main 
results achieved are the following: 

1. Although there exists no origin, where R would be 
zero permanently, a set where M = can exist. At 
this location, R, t is constant (section [TTTJ) . 

2. The whole spacetime is both future- and past- glob- 
ally trapped (section IV). 
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3. The geometrical interpretation of the (x, y) coor- 
dinates in a constant- (t, z) surface was clarified in 
Sec. El Contrary to an earlier claim [l[, this sur- 
face consists of just one sheet, doubly covered by 
the (x, y) map. 

4. The geometries of the following surfaces for the 
metric (|5.1[) were shown in illustrations, all of them 
in Sec. |VD 



(a) z = constant, <p — for (|5.1|) in Figs. H]to|31 

(b) The collection of R(t, z) curves in Fig. [5] 

(c) t = constant, tp = for (|5.1[) in Figs. [6] and [8] 

It turned out that the surfaces listed under (a) are 
locally isometric to ordinary surfaces of revolution 
in the Euclidean space (in special cases to a plane 
and a cone) when E > 1, but cover the latter an 
infinite number of times. When 1/2 < E < 1, 
they cannot be embedded in a Euclidean space even 
locally. The values E < 1/2 are prohibited by the 
spacetime signature. 

The time evolution of the surfaces under 3(c) was 
illustrated in Figs. 151 and ITU1 

5. For the general metric (I2.7P - (|2.8I) . the geometry of 
the surfaces t = constant, ip — was investigated in 
Sec. IVIII and shown in Fig. [TTJ The other surfaces 
listed above are the same as in the hyperbolically 
symmetric case (|5.1[) . 



6. In Sees. I Villi to IXII a detailed analysis was car- 
ried out of the relation between the mass function 
M(z) and the sum of rest masses in a volume A4(z) 
in the quasi-spherical Szekeres model. The func- 
tion M(z) represents the active gravitational mass 
within a sphere of coordinate radius z, while M(z) 
is the sum of rest masses of particles contained in 
the same volume: 



M 



\93\pd3X, 



(14.1) 



mass distribution cause difficulties that were not 
fully resolved. It has only been demonstrated that 
the average value of M/V over the whole space 
t = to is determined by the average value of M/Vb, 
where Vq is the flat space limit of V. 

This problem requires further investigation, but it 
is hoped that the results achieved here will be of 
use for that purpose. 



Appendix A: The curvature tensor for the metric 



The formulae given below are the tetrad components 
of the curvature tensor for the metric (|5.1[) . The tetrad 
is the orthonormal one given by 



df, 



R. 



V2E-1 



dz. 



c 2 = Rdti, e 3 = Rsinhddip, 



(Al) 



with the labeling of coordinates (t, r, <p) — 
(x° , x 1 , x 2 , x 3 ). The components given below are 
scalars, so any scalar polynomial in curvature compo- 
nents will be fully determined by them. Since they do 
not depend on they have no singularity caused by any 
special value of □ 



Ro 



2M M, 



R 3 R 2 F' 
1 



-R(J202 — ^?0303 — 2 -^2323 - 

M M, 
R\2\2 — R\Z\Z — TP? 



R? R 2 F' 



(A2) 

■w (A3) 

(A4) 



The formulae in both appendices were calculated by 
the algebraic program Ortocartan [29|, [3(| ■ 



where V is any volume in a space of constant t — to, 
53 is the determinant of the metric in that space 
and p is the mass density at t = to- The relation 
(|8.3[) follows in the limit when V is the volume of 
the whole space t — to. The calculations in Sees. 
IVIIII to IXII demonstrated how to calculate (|14.1|) 
within various relevant volumes. 



7. In Sees. IXIII and IXIIII calculations analogous to 
those from Sees. IVIIII - IXII were carried out for 
the quasi-hyperbolic Szekeres models. The aim 
was to interpret the function M{z) in this case by 
identifying the volume in which the active gravi- 
tational mass is contained. Integrals analogous to 
(114.11) can be calculated, but the full analogy with 
the quasi-spherical models follows only in the (hy- 
perbolically) symmetric case. In the general case, 
terms arising from the dipole component of the 



Appendix B: The curvature tensor for the metric 



The formulae given below are the tetrad components 
of the curvature tensor for the metric (|2.8p with e = — 1. 
The tetrad is the orthonormal one given by 



c° = dt, e 1 = 



F 



dz, 



V2E-1 

2 R , 3 R , 
c = — da;, e° = — dy, 

with the labeling of coordinates (t, z, x, y) 
(a; ,^ 1 ,^ 2 ,^ 3 ), where £ is given by (|2.7p and 



(Bl) 



def 



F — R, z —R£, z / £. 



(B2) 
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The components given below are scalars, so any scalar 
polynomial in curvature components will be fully deter- 
mined by them. 



-Roioi = 
-R0202 = 
-R1212 = -R1313 



2M 3M£, 2 



M,. 



R 3 

F0303 



R 2 £F R 2 F' 
1 M 

^1-2323 " 



2 
M 



i? 3 ' 
ZM£, Z M,, 



R 2 £F R 2 F' 



(B3) 
(B4) 
(B5) 



Note that these reproduce (|A2|) - ([A4]) when £ , z = 0. 

We wish to find out whether the sets £ = and £ — > 00 
are curvature singularities. For this purpose it is useful 
to introduce the coordinates <p) by (|5.3p . Since the 
quantities (|B3[) - (IB5|) are scalars, we only need to sub- 
stitute (|5 .3[) in them. The two suspected sets become 
i? — > 00 and $ = 0, respectively. After the transforma- 
tion we have 



2sinh^(^/2)' 



(B6) 



£ )2 



2sinli (i?/2) 



[1 - 2cosh 2 (z9/2)] 



- coth(tf/2) (P,z cos <p + Q, z sin 95) . (B7) 



The only quantity in (|B3[) - (IB5|) that depends on $ is 
£, z Using dHSJ) - (E3) we easily find 



,. £,z 1 

lim — — = — . 

■d^co£F R' 

lim — ^ = 



-FS, Z /S 



(B8) 
(B9) 



The loci where these can become infinite do not depend 
on ??. Hence, i? — > 00 and = are not curvature singu- 
larities, and neither are £ = or £ ^ 00. □ 
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